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Abstract 



Final states involving tau leptons are important components of searches for new particles at the Large 
Hadron Collider (LHC). A proper treatment of tau spin effects in the Monte Carlo (MC) simulations 
is important for understanding the detector acceptance as well as for the measurements of tau polar- 
ization and tau spin correlations. In this note we present a TauSpinner package designed to simulate 
the spin effects. It relies on the availability of the four- momenta of the taus and their decay prod- 
ucts in the analyzed data. The flavor and the four-momentum of the boson decaying to the T~ T + or 
T ± V pair need to be known. In the Z/y* case the initial state quark configuration is attributed from 
the intermediate boson kinematics, and the parton distribution functions (PDF's). TauSpinner is 
the first algorithm suitable for emulation of tau spin effects in tau-embedded samples. It is also the 
first tool that offers the user the flexibility to simulate a desired spin effect at the analysis level. An 
algorithm to attribute tau helicity states to a previously generated sample is also provided. 
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1 Introduction 



Tau leptons are an excellent signature with which to probe new physics at the LHC. As the heaviest 
leptons, they have the largest coupling to the Higgs boson both in the Standard Model (SM) and the 
Minimal Supersymmetric Standard Model (MSSM). Their short-enough lifetime and parity-violating 
decays allow for their spin information to be preserved in the decay product kinematics recorded in the 
detector. They are therefore the only leptons suitable for measuring longitudinal polarization and its 
correlations, which provide important constraints on the nature of the observed resonance. 

In this note we present a TauSpinner package, which is a MC program designed to generate tau spin 
effects in any tau sample provided their origin is known. This tool has two important applications: 

Data driven analysis of Z/y* ->■ t T + and W ± — > T ± v 

backgrounds. In the case of the algorithms for embedding taus on measured light lepton sam- 
ples HI the TauSpinner may represent a third step. The first step is to construct the tau four- 
momenta from the four-momenta of the measured lighter leptons (and accompanying photons) 
while the second step comprises the decay of unpolarized taus usin^ 

e.g. Tauola (2l[3l|4l. Omission of the last step leaves the kinematics of the tau decay particle 
different from what is expected from the polarized taus that appear in nature. This could lead to a 
mis-modeling of the shapes of various observables and hence a mis-measurement of the acceptance 
of a given set of cuts. 

MC studies. Since the tau spin effects affect the overall acceptance of the taus in the detector, a proper 
treatment of these effects is of great importance for the interpretation of the results as well as 
feasibility studies of new models. It is also key for measurements of tau polarization and tau 
spin correlations. The TauSpinner algorithm allows one to create samples with different tau 
polarizations from an initial sample by re-weighting events to give the desired distributions as a 
function of the decay mode of the tau. Furthermore, the helicity states of the taus can be attributed 
to the previously generated sample (with spin effect included or introduced with the TauSpinner 
weight). 

This letter is organized as follows. In Sec. [2] we present the algorithm while in Sec. [3] we discuss its 
performance as compared to the standard Tauola MC package. The results are summarized in Sec. [4] 
All relevant technical details are provided in App. A and B. 

2 The TauSpinner algorithm 



The TauSpinner algorithm relies on a leading order approximation in which spin amplitudes are used 
to calculate the spin density matrices for hard 2 — > 2 or 1 — > 2 Born level processes ||5l |U 13 HI 0. 
TauSpinner is constrained to the longitudinal tau spin degrees only. It starts with identifying the flavor 
of the intermediate boson: W ± , Z/y*, H or H ± . The information on the four-momenta of the outgoing 
taus and their decay products as well as the intermediate boson four-momentum is then used to determine 
the polarimetric vectors. The longitudinal tau polarization (P t ) is randomly generated as specified in 
Tab. [T] and set to ±1 which correspond to pure tau helicity states. Probability of the helicity states is a 
constant for taus originating from the W ± , H ± and H bosons. In Z/y* — > T~T + events, the probability, 
denoted , is a function of the T~ scattering angle, 6, and the center of mass squared of the hard process, 
s. At the Born level, and in the ultrarelativistic limit, is, following notations of [5], given by: 

'in order to enhance statistics for the spin analysis, the decay of unpolarized taus can be performed multiple times. 
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da 



a (5,COS0,P T = 1) 
acos8 v 1 I / 

d^se(^cos0A = l) + a^e(^cos0,F T ,= -l 



where 



(s, cos , P t ) = ( 1 + cos 2 )F (j) + 2 cos 0Fi (j) - P % [( 1 + cos 2 )F 2 (5) + 2 cos 6F 3 {s)}, (2) 

d cos \j 

and fj (5) 's are four form factors which depend on the initial and the final state fermion couplings to the 
Z boson and the propagator. The dependence on the T~ longitudinal polarization P x is taken into account 
(note that in this case P % = ±1, also P t] = P Zl = P T ). In [5] the probability is calculated using the 
information on the initial state quarks stored at the generation level. In TauSpinner the initial state 
quark configuration is attributed stochastically from the intermediate boson kinematics and PDFs in the 
following steps: 

1. The invariant mass of the Z/y* is calculated from the intermediate boson four- vector. Note that it 
does not need to coincide with the sum of T~ T + four-momenta as it may include photons of final 
state bremsstrahlung. 

2. The scattering angle, cos 0, is calculated in the t~t + 

pair rest frame from the angle between the direction of the first beam ( 1 ,0,0, 1 ) boosted to this frame 
and the direction of the T + or from the angle between the direction of the second beam (1,0,0,-1) 
boosted to this frame and the direction of the T~ . In the final step the average of the two 

3. The fraction of momenta taken by partons of the proton: x\ and X2 are resolved from constraints 
x\X2Eq M = s and (jci — X2)Eqm = Pz, where £"cm 

and p z denote the collision center of mass energy and the longitudinal component of the Z/y* mo- 
mentum, respectively. 

4. Probabilistic choice on the basis on the leading order 2—^2 Born level cross sections and PDF's 
is performed to attribute the flavors to the incoming quarks and sign of 0. 

The probability is calculated as a weighted average over all possible initial state quark configurations. 
This solution is only implemented for the case of proton-proton (pp) collisions. For other types of 
collision events the functionality of TauSpinner for Z/y* — > T~T + process is restricted to longitudinal 
spin correlation only (p^ =0.5). 

Note that the calculation of the momentum fractions x\ and X2 is performed within a collinear approxi- 
mation where the initial state radiation quarks and gluons are assumed to have no transverse momentum. 
Furthermore the final state radiation photons are taken into account at collinear level only: they are omit- 
ted from the boson decay vertex. Any deviation from the energy-momentum conservation in the boson 
decay vertex is attributed to the presence of the photons. As shown in Sec.[3]these approximations have 
practically no impact on the implemented spin effects. 

2.1 Calculation of the spin weight 

The outcome of running the TauSpinner program is a spin weight attributed to each event separately. 
In the T V final state, for any decay of a polarized tau, the spin weight is defined as [2j|3l|4l: 
2 In (Tp 1 this angle is referred to as 9* is taken 1101 
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Table 1: Probability for the configuration of the longitudinal polarization of taus from different origins |5|. 



Origin 






Probability 


Neutral Higgs bosons: H 


+ 1 


- 1 


0.5 




-1 


+ 1 


0.5 


Neutral vector boson: Z/y* 


+ 1 


+ 1 


P? 


-1 


-1 




Charged Higgs: H^ 1 


+ 1 




1.0 


Charged vector boson: W ± 


- 1 




1.0 



w T = 1 + s-h (3) 

where s is the tau polarization vector, and h is the polarimetric vector constructed using hadronic currents. 
In our W ± and H ± decays the exact expression reduces to: 

wt = 1 + sign h z (4) 

where sign equals one for left-handed taus from W ± bosons and minus one for right-handed taus from 
the charged Higgs boson. h z is the z component of the polarimetric vector. 

In the t~t + final state, the weight is defined as B21 EJ SI : 

w T = Rqfttf (5) 

where /?y is a matrix describing the full spin correlation between the two taus as well as the individ- 
ual spin states of the taus. h 1 and ti ] are the time (?) and space (x, y, z) components of the two taus' 
polarimetric vectors. The t component of h and R tt are, by convention set to 1. 

Neglecting the transverse spin degree, and in the ultrarelativistic limit, for Z/y* — >■ T~T + events the 
expression reduces to: 

wj = 1 + sign h z +h z - + P T h z + +P z h z - (6) 

where P T denotes the polarization of the single tau in a mixed quantum state. Within this approximation, 
P x is a linear function of the probability : 

P* = 2p?-1. (7) 

In an event of a neutral and spin zero Higgs boson decaying to t~t + , the expression simplifies to: 

wt = 1 + sign h z +h z - (8) 

The sign equals one for the Z/y* boson and minus one for the neutral Higgs boson, reflecting the opposite 
spin correlations in the two samples. 

Each tau decay channel requires a distinct method to calculate the polarimetric vector [6 ]. The tau decay 
modes implemented in TauSpinner are listed in Tab. [2] For the remaining channels, involving five pions 
in the final state and multi-prong decays with kaons, the effect of tau polarization is neglected and h z ± is 
set to zero. 

2.2 Application of the spin weight 

The event weight can be used at the analysis level for: 

Simulating tau spin effects e.g. in a sample generated without spin effects. The event weight equals 
wt- It takes values between (0,2), except for the case of Z/y* — > T~T + when the range is (0,4). 
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Table 2: Summary of tau decay modes implemented in the TauSpinner. Branching fraction is given for each 
decay mode ifTTIl . 



Tau decay mode 


Branching fraction % 


e v e v T 


17 R5 




17.36 




10.91 




25.51 




9.29, 9.03 (incl. co) 


K v 


0.70 


K-7i°V, n-K°v 


0.43, 0.84 


n~n + K~n°v 


4.54 (incl. co) 


7l~K K K°V 


1.04 


Other 


2.5 



Removing tau spin effects from a sample generated with spin effects. The event weight equals 1 /wj. 
Then, it is greater than zero with no upper limit. 

Reverting tau spin effects in a sample generated with certain longitudinal tau polarization (and/or cor- 
relations) to the different one. The weight equals (2-wj)/wt for the z ± v final state corresponding 
to the W ± ^H ± or H ± ^W ± replacement, and w T (H)/w T (Z/f) or w T (Z/f)/w T (H) for the t^t+ 
final state corresponding to Z/y*—)-// or H—t-Z/y* replacement. It is greater than zero without an 
upper limit. 

2.3 Attributing tau helicity states in the Z/y* — > t T + decays 

The helicity states of the taus are attributed stochastically by comparing a random number with the 
probability, p^ wj{p^= 1)/wt, of the right-handed configuration to occur. This method is valid for 
events generated with spin effects or emulated using the spin weight provided by the TauSpinner. An 
average over all possible initial state configuration is taken. Presentation of the method principle can be 
found in Ref. 0. 

3 Performance of the TauSpinner algorithm 

The performance of the TauSpinner algorithm is studied based on the MC simulated Z/y* — > T~T + and 
— > T ± v events. The samples were generated using the general purpose event generator Pythia lfT2l 
assuming the pp collision at center of mass energy of 7 TeV. The taus were then made to decay using the 
Tauola package. Two sets of events were simulated: 

No spin effects. In these events the taus were decayed by Tauola as if they had been produced with no 
polarization. They were used for emulation of spin effects later, with the TauSpinner package. 

Tauola. In these events tau spin effects were properly accounted for at the time of tau event generation 
and decay. 

Event generation and simulation of the spin effects in the TauSpinner were performed with the modified 
LO partem distribution function (PDF) MRSTLO* Q3). 

3.1 Simulation of the tau polarization in t ± v final state 

For all tau decay modes, the main observables that are sensitive to the tau polarization are the tau mo- 
mentum fraction taken by the hadronic system lfl4l . x, and the relative difference between the charge and 
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(a) Fraction of the tau momentum taken by the hadron in (b) Relative difference between the charged and the neu- 
the combined K~ v and K~ v channels. tral energy in the K~ 7t l) v channel. 



• No spin effects 

* W*-*^ (TauSpinner) — 
A hF-»T-v (TauSpinner) 

: W*-^ (Tauof 




B I'D 

1 

0.8 
0.6 
0.4 
0.2 




• No spin effects 

+ W 1 -*!^ (TauSpinner) 

ir hP-^v (TauSpinner) 

I W*-^ (Tauola) 




(c) Relative difference between the charged and (d) Relative difference between the charged and the neu- 
the neutral energy in the combined jr _ jr JT v and tral energy in the 7C~7i°v channel. Decays of taus in the 
7t~7i~n + v channels. "no spin effects" sample were generated using Pythia 

6.425. 



Fig. 1: Comparison of spin effects simulated with Tauola and TauSpinner for various observables. 



the neutral energy in the tau decay, T JT51I . 

Plots a)-c) in Figure [T] demonstrates the performance of TauSpinner for the channels where the tau 
polarization is extremal. The observable x is plotted in Fig. 1(a) for the combined 7i~v and K v chan- 

for the 7i~7i°v and combined n~K°K Q v and 



nels. The observable T is plotted in Fig. 1(b) 



and 



1(c) 



7i~7t~7t + v channels, respectively. The sample with no spin effects refers to W 
events generated with flat P T value. The H — > T V and W ± — > z ± v configurations were obtained by ap- 
plying an appropriate spin weight to the sample with no spin effects. The weighted observables exhibit 
the expected behavior indicating a proper implementation of the P T in the TauSpinner package. 

Note that the TauSpinner algorithm requires that the decays of unpolarized taus via p or a\ preserve the 
spin correlation between the production and decay of these mesons. This feature is missing if the taus 
are made to decay using Pythia 6.425. Although a bulk of spin effects are reconstructed, a significant 
systematic uncertainty arises as demonstrated in Fig.[T]d). 



3.2 Simulation of tau polarization in Z/y* — )> z T + events 

InZ/7* — > t~t + events, the single tau polarization P t depends on the intermediate boson virtuality, the 
flavor of the incoming quark and the scattering angle. The complexity of this dependence is shown in 
Fig.|2]where the tau polarization is drawn as a function of cos 8, for two different ranges of the invariant 
mass of the boson. In these plots, the quark flavor configuration is fixed at the level of simulation of the 
spin weight in the TauSpinner. 

In the next step, a degree of polarization is studied for different configurations of the incoming quarks. 
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(a) The up quark aligned along the positive z axis. (b) The down quark aligned along the positive z axis. 



Fig. 2: Tau polarization as a function of cos 6. 




(a) The up and down quarks aligned along the positive z (b) The up and down quarks aligned along the negative z 
axis. axis. 

Fig. 3: Fraction of the tau momentum taken by the hadron in the K~v and K~v channel in Z/y* — > t~t/ + events. 
Taus with negative charge emitted in the forward hemisphere is chosen. 



The virtuality of the intermediate boson is constrained to lie within ±3 GeV of the Z boson mass and both 
taus are set to decay to %~v or K v. The forward-backward spin asymmetry is accessed by choosing the 
T~ to be emitted in the forward region by requiring the longitudinal momentum of T~ to be greater than 
that of T + . Figure |3] shows the observable x for the up and down type quarks entering along the positive 
or negative z axis. The results are consistent with those of reference [5], and therefore reassure a proper 
transmission of spin effects from the hard process to the tau decay products. 

In the last step, the degree of polarization is studied inclusively for all initial state quark configurations 
and the results are compared to those simulated using the Tauola package. The Z/y* virtuality is con- 
strained to lie within ±3 GeV of the Z boson mass. The fraction of the tau momentum taken by the 
hadron in the 7i~v channel is plotted in Fig. |4j Expecting the Z boson to be emitted in the direction 
of the incoming quark (as compared to the direction of the anti-quark from the sea) these plots can be 
compared to those in Fig. [3] 

Figure [5] shows the fraction of the tau momentum taken by the hadron in the combined tt~ v and 
K v channels for the left-handed and the right-handed taus from the Z/y* — > T~T + decays. Helicity 
states are attributed by the TauSp inner in the sample without and with spin effects simulated at the gen- 
eration level. The distributions exhibit the expected shapes for properly attributed pure helicity states. 
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Z— in* (Tauola) 
Z-hV (TauSpinner) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



(a) Z boson emitted in the forward direction. 



(b) Z boson emitted in the backward direction. 



Fig. 4: Fraction of the tau momentum taken by the hadron in the % v channel. Taus with negative charge emitted 
in the forward hemisphere are chosen. 




X 



Fig. 5: Fraction of the tau momentum taken by the hadron in the combined % v and K v channels. Xr and i£ 
denote right-handed and left-handed taus, respectively. 
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(a) The z s variable described in text in the channel where (b) Visible mass of the two hadrons in the channel where 
both taus decay to lt~ v or K~ v. both taus decay to n~ v or K~ v. 



Fig. 6: Tau spin correlation observables. 



3.3 Validation of spin correlation in the z T + flnal state 

Two observables were demonstrated in to be appropriate for studying tau spin correlations in the 
T~T + final state: the invariant mass of the hadronic system, m v i s , and the z s variable. The latter is 
a signed surface in the x + - x~ plane, between lines: x + =x~ and x + =x~+a. x + and x~ denote the 
fraction of tau momenta taken by the hadrons. The two observables are plotted in Fig. [6] in the channel 
where both taus decay to 7i~v or K v. The sample with no spin effects refers to Z/y* — > t~t + events 
generated with =0.5 and no spin correlation. The proper H — > %~ z + and Z/y* — > %~ T + configurations 
were obtained by applying an appropriate spin weight to this sample without spin effects. The weighted 
observables exhibit the expected behavior reassuring a proper implementation of the spin correlations in 
the TauSpinner package. 



4 Summary and outlook 

The TauSpinner package designed to emulate tau spin effects has been introduced. The algorithm 
is limited to the leading order accuracy and the longitudinal spin degrees only. As compared to the 
algorithm developed in Ref. [5], the functionality is extended to allow one to estimate the spin effects 
when the information on the incoming quarks entering the hard process is not available. This is therefore 
the first algorithm suitable for emulating the spin effects in the tau-embedded samples. The intrinsic tau 
polarization arising from parity violation in the weak interactions is attributed on the basis of intermediate 
boson kinematics and PDF's. For the Z/y* — > t~t + process its complete functionality is limited to pp 
collisions and otherwise restricted to conservation of total angular momentum. 

Comparisons of Tauola and TauSpinner for various spin observables presented in Sec. [^demonstrate 
that no additional systematic uncertainty in simulation of tau spin effects has been introduced in the 
TauSpinner. This claim should, however, be validated in a sample where the number of high transverse 
momentum jets is enhanced. 

A complete discussion of the theoretical uncertainties is common for the Tauola, the TauSpinner 
and, to a large extent the KORALZ [ 8 ] MC programs. Some aspects of these uncertainties have been 
addressed through the work published in Refs. EJ [3j III El HH . A rigorous evaluation of theoretical 
effects and in particular a comparison with results based on the next-to-leading order matrix element 
calculations simultaneous for scattering processes and for spin density matrices is referred to a future 
work. 

Although the current version of the TauSpinner is restricted to the longitudinal spin degrees, the frame- 
work is prepared to simulate the complete spin effects. This extension is planned for the future version 
of the algorithm. The code of TauSpinner is publicly available, with all relevant details given in App. 
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A and B. 
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A Requirements for data files 



Data files generated by a MC event generator (or constructed using the tau-embedding method) need to 
fulfill the following requirements: 

1 . Four-momenta of the intermediate boson, the taus and the flavor and the four- momenta of the tau 
decay products need to be available. 

2. Flavor of the intermediate boson needs to be available or set by the user. 

3. For all types of hard processes, but Z/y* — > T~T + , the four-momenta can be defined in an arbitrary 
but common frame. For the Z/y* — > z~z + process, the four-momenta have to be given in the 
laboratory frame in order to be consistent with the PDFs. 

4. The four-momenta of the taus and their decay products need to be known with sufficient precision 
in order to ensure numerical stability of the algorithm. Six significant digits are recommended. 

Note that different MC generators may store the truth information in different ways. It is the responsi- 
bility of a user to make sure that all these requirements are fulfilled. 

B Public version 

A generic version of the package can be found in iPTTI . The main code is written in C++ and relies 
upon two libraries: Tauola and LHAPDF [18]. A method for reading input information stored using the 
HepMC [ 19] format is prepared. Support for any other input format is available upon request. 

B.l Project organization 

The TauSpinner package is organized in the following manner: 

- src/tau_reweight_lib . c , src/tau_reweight_lib .h - the core of the algorithm. 

- src/Tauola_wrapper . h - wrapper for TAUOLA FORTRAN routines. 

- src/SimpleParticle .h - definition of class SimpleParticle used as a bridge between the 
event record (or data file) and the algorithm. 

- src/Particle .h - definition of class Particle used for boosting and rotation of the particles. 

- src/read_particles_f rom_TAU0LA . c , src/read_particles_f romJTAUOLA . h - interface to 
the HepMC : : IO_GenEvent data files used by the example program. 

- README - a short manual. 
B.2 The algorithm sequence 

The TauSpinner takes the following sequence of steps: 

Initialization of Tauola. It is ensured by invoking: 
Tauola: : initialize () ; 

Initialization of TauSpinner. It is performed by executing: 

void initialize_spinner (bool Ipp, int Ipol, double CMSENE) 

where the argument Ipp passes the information on the type of collision events (Ipp = true sets 
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pp collisions), Ipol passes the information on the spin effects included in the input sample (Ipol= 
0, 1, 2 corresponds to no spin effects, complete spin effects and spin correlations only, respectively) 
and CMSENE sets the collision center of mass energy. 

Reading the data files. Information on the four-momenta and the flavor of the boson, the final state taus 

or tau and neutrino pair and the tau decay products is filled and stored in instances of SimpleParticle 
class by the use of the function: 

void readParticlesFromTAUOLA_HepMC(HepMC: : ICLGenEvent &input_f ile, SimpleParticle 
feboson, 

SimpleParticle &tau, SimpleParticle &tau2, vector<SimpleParticle> &tau_daughters , 
vector<SimpleParticle> &tau2_daughters) . 

This function should be modified if the input files are not in the HepMC : : ICLGenEvent format. 

Calculation of the spin weight. It is performed by the use of the following functions: 

double calculateWeightFromParticlesWorHpn (SimpleParticle feboson, SimpleParticle 
&tau, 

SimpleParticle &tau2, vector<SimpleParticle> &tau_daughters) forthelV 111 — > TV and 
H — > t~t + processes 

double calculateWeightFromParticlesH (SimpleParticle feboson, SimpleParticle &tau, 
SimpleParticle &tau2, vector<SimpleParticle> &tau_daughters , 
vector<SimpleParticle> &tau2_daughters) for the Z/y* — > T~T + and H — > T~T + pro- 
cesses. 

Attributing tau helicity states. For Z/y* — > T~T + process, the tau helicity states are attributed at the 
stage of calculation of the spin weight. The information can obtained by calling getTauSpinO 
function. 

B.3 Calculation of the spin weight 

For the T^v final states, the spin weight is calculated in the following steps: 

1. The parent boson, the tau, the tau neutrino and the list of tau daughters are identified and boosted 
to the T^v rest frame in which the tau is aligned along the z axis. 

2. The tau daughters are boosted to the tau rest frame. Two angles of spacial orientation of the 
neutrino from the tau decay, theta2 and phi2, are calculated and stored. The tau daughters are 
rotated by these angles in order to align the neutrino along the z axis. 

3. The Tauola decay channel is identified. 

4. The Tauola FORTRAN subroutine is called to perform calculation of the polarimetric vector h. 

5. The polarimetric vector h is rotated back using the theta2 and phi2 angles. 

6. The spin weight is calculated using eq.|4]and returned to the main program. 

For the T~T + final states, the spin weight is calculated in the following steps: 

1 . The parent boson, the taus and their tau daughters are identified and boosted to the T~ T + rest frame 
in which the taus are aligned along the z axis. 
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2. For each tau: 

Its identified daughters are boosted to its rest frame. Two angles of spacial orientation of the 
neutrino from the tau decay, theta2 and phi2, are calculated and stored. The tau daughters are 
rotated by these angles to align the neutrino along the z axis. 

The Tauola decay channel is identified. 

The Tauola FORTRAN subroutine is called to perform calculation of the polarimetric vector 

h. 

The polarimetric vector h is rotated back using the theta2 and phi2 angles 

3. In case of the Z/y* — > t~t + decays: 

The probability is calculated using eqs[j][2] 

The spin weight is calculated using eq.[6]and returned to the main program. 

4. In case of the H — > t~t + decays: 

The spin weight is calculated using eq.[8]and returned to the main program. 

B.4 The LHAPDF library wrapper 

The evolution of the PDF's is invoked from the wrapper for PDF's: 
double f (double x, int ID, double SS, double cmsene) 

where function f calls the evolution function xfx(x, SS, ID) lfl8Tl . The PDF sets need to be avail- 
able locally. They can be obtained from the LHAPDF project website. 
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